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The calculation of quantum defects in an isoelectronic 


series based on polarization concepts 


By Howarp A. Rosprinson 


It is often of importance, when analyzing a spectrum, to be able to predict 
a given quantum defect in order to obtain at least approximate absolute term 
values. The classical method for obtaining such absolute values is of course 
to fit terms of running total quantum number » to a Rydberg or a Ritz 
series. Very often, however, a sufficient number of such terms cannot be found 
and particularly in cases where large values of n and high stages of ioniza- 
tion are involved, experimental accuracy of wave length measurement may not 
be high enough to give particularly accurate term values. It would, therefore, 
be convenient to have methods for extrapolating effective quantum numbers 
along an isoelectronic sequence. A method for accomplishing this in the case 
of non-penetrating orbits is given in this paper. 

In order to make as much use as possible of methods which require the 
least amount of computation one naturally looks for means by which the least 
number of significant figures may be used to get the highest degree of ac- 
curacy. Thus one is tempted to consider first not the absolute term values 
themselves but rather the deviation of the absolute values from the cor- 


responding hydrogenic values. Thus if 7y = — a represents a hydrogen 
like term 
R2 Rlz + (x — o)f* 
bust ne pote he 


will represent a non hydrogen like term if z—1= the stage of ionization, 
Z=2z+x2-=the atomic number and o =the screening constant. One can then 
write 


R 
1 — Tq =—2(x—0) ay — (x0) p= Art B 
Nn nv 
where 
R »R 
A=—2(x—0)-4, B=—(%—a)* 


Now unfortunately B and A are not true constants in most isoelectronic se- 
quences and thus 7’ — Ty is not strictly linear. This relation thus does not 
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lend itself too readily to extrapolation use inasmuch as it is actually as simple | 


to estimate changes in * with z as it is to estimate changes in o with z. 
A relationship between n* and o is easily derived. Since 


ik (@t+x—o ~— 2 
po en ee 
the equation 
Z Ne 
2 ae) 


follows immediately. It is of interest to examine this relationship for possible 
minima. In the usual manner it can be shown that n* will be a minimum 
when 


d(% — o) doy nile Or 


dz a ok ITZ 


It has been known for many years that there are in general two types of 
curve found when n* is plotted against z. In the case of penetrating orbits, 
which are usual with s and p electrons, no minimum is found and the curve 
increases continuously with z. In the case of non-penetrating orbits, however, 
the curve goes through a clear minimum. Examples of these two types of 
behavior are given in Table I and Table II.? 

It is clear that in those cases where one has a minimum ahead! 


Az 


Thus the presence of a minimum is caused by the value of o. 


is equal 


toes 


This in turn is the justification for the concept of penetrating and non- 
penetrating orbits. 

It is possible by means of another method of attack to obtain extrapola- 
tion functions of greater utility particularly in the case of non-penetrating 
orbits. Previous work? has shown that quantum defects in non-penetrating 
orbits are largely due to polarization effects. Simple theory developed in the 
cited references leads to the result that, due to self polarization within the 
atom, a better approximation to the true term energy is given by 


Ww R2 ee (5) 


he n> Qhe\r4 


* Data in this paper were taken from the following sources : 


The ions of Li through F: From various papers by Epimn as given in Atomic Energy 
Levels, National Bureau of Standards Circular 467, Washington D.C. 

The ions of Na, Mg, Al: SépERQvist: Arkiv for Matem. Astr. och Fysik 30 A, No 11 
(1944), 32 A, No 19 (1946). 

The ions of Si and §: Ferner: Arkiv for Matem. Astr. och Fysik 36 A, No 1 (1948). 

The ions of P: Unpublished Data of the Author, 


2 PauLina: Proc. Roy. Soc. A 114, 181, 1927. — Mayer and Mayer: Phys. Rev. 43, 
605, 1933. — Van Vurck and WuirELaAw: Phys. Rev. 44, 551, 1933. —- Paurine and Goup- 
smit: The Structure of Line Spectra, McGraw-Hill Book Co. N.Y. 1930. 
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O I isoelectronic sequence 2 p*%P, (penetrating orbit). 


Table I 
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Ton ne (% — 0) MOY) o lb 
Az Z 
ORL 1.0005 0.9990 6.0010 0.9990 
0.2084 : 
thea Ge 1.2471 1.2074 5.7926 .6037 
(1092) 
Ne UT (1.390) (1.3166) (5.6834) -4389 
(0754) 
Na IV 1.4837 1.3920 5.6080 3480 
.0526 
Mg V 1.5517 1.4446 5.5554 2889 
-0386 
Al VI 1.6036 1.4832 5.5168 .2472 
0295 
Si VII 1.6446 1.5127 5.4873 2161 
.0159 
Pe Vir 1.6774 1.5386 5.4614 1923 
.0423* 
S Ix 1.7051 1.5809 5.4191 1757 
* The increase here is possibly due to experimental inaccuracy. 
Table IT 
C I isoelectronic sequence 3d°D, (non-penetrating orbit). 
Ton n* (% — 0) Aso) o ae 
Az z 
Oh AL 2.9556 0.01502 4.98498 0.01502 
0.03744 
N II 2.9236 -05246 4.94760 .02623 
-02998 
O Il 2.9198 .08244 4.91756 .02748 
.02844 
ne LY: 2.9191 11088 4.88912 02772 
(.017) 
Ne V (2.925) .128 (4.872) (.0256) 
(.015) 
Na VI 2.930 14334 4.85666 02839 
01416 
Mg VIL 2.934 .15750 4.84250 02250 
.01130 
Al VIII 2.938 .16880 4.83120 02110 
00859 
Si IX 2.942 T7389 4.82261 01971 
.01281* 
1p 2 2.944 .19020 4.80980 01902 


* The increase here is possibly due to experimental inaccuracy. 
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Mee ners 1 if 
where a is the polarizability and (;:) the average value of the reciprocal fourth | 
power of the radius of the electronic orbit, is given by 


aj n® (1 + 3/,) (1+ 1) (1 + 3) (2) 74) 


Since for extrapolation purposes one is primarily interested in the energy varia- | 
tion with Z (or z) the question arises as to how a varies with Z. PAauLine 
shows that- a= Z—s) and remarks that the polarization screening constant © 
S is apparently independent of Z in the Li I isoelectronic sequence. There — 
are, however, other smaller effects which come into play, such as interaction _ 
due to some penetration and interaction due to quadrupole effects. Further- 


; ae 1 
more published values of a show more nearly a variation with Zs Such 
effects must also be taken into consideration. One may write 

W R2* ez =e | ice es 
he = =(n—A) no 2he Ve) fie 
or, substituting 
fete 
ane 
W ze Ge ae ii 
2 ay (n — A)” 2 ay n* 2 \r4 
Since, assuming 4 <n, 
i. ae ( pene 
(n—Al® wr? n ne 
3 
1 4 
A= Sr (") = K' az? + higher terms 
22° \r 


where 
K’ = 1.1245 X 10%* for 2> electrons 


= 04998 x 1074 for 3d electrons. 


It is logical in view of the above mentioned dependence of a on z to try an 
extrapolation function of the form 


a? 
~ f(z) 
64 


Table III 
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Li I isoelectronic sequence. Term 3d?Ds), 
(z + 0.527)® + 2.63402 + 0.9413 = 0.10133 22/A 


Spin orbit correction = 0.0191 24 


n* cale. 


2.998580 
2.998186 
2.998270 
2.998445 
2.998615 
2.998763 
2.998886 
2.998990 
2.999077 
2.999151 
2.999215 
2.999270 
2.999318 


Term eale. 
em? 


12 203.62 
48 428.22 
109 858.3 
195 282.0 
305 095.6 
439 296.6 
597 885.2 
780 858.4 
988 220.2 
1 219 966 
1 476 099 
1 756 618 
2 061 524 


Spin orbit 


Total term 
cale. 


em? 


12 203.64 
48 828.53 
109 859.8 
195 286.9 
305 107.5 
439 321.4 
597 931.1 
780 936.6 
988 345.5 
1 220 157 
1 476 379 
1 757 014 
2 062 070 


Term obs. 
em + 


12 203.64 
48 828.5 
109 859.9 
195 286.9 
305 105.9 
439 322.7 
597 934.0 
988 316* 
1219 930* 
1 476 360* 
1 756 930* 


(Obs.—cale.) 
= 


em 


ok 
i= 


* These terms were found absolutely by extrapolation of n* not by calculation of a Ritz 


series. 
Table IV 
Li I isoelectronic series. Term 2 p?P»), 
z* + 377.681 z? + 929.960 2? + 647.511 z + 441.808 = 97.513 2?/A 
Spin orbit correction = 0.091314 (z + 0.1). 
* Term calc. Spero acco) term Term obs. | (Obs.—cale.) 
Zz n* cale. 4 corr. eale. aA aa 
em ai =| em em 
em cm 

Li I 1.959335 28 582.65 0.13 28 582.78 28 582.76 — (0.02 
Be IT 1.954080 114 948.4 1.8 114 950.2 114 946.3 = es 
B Ii 1.958273 257 530.1 8.4 257 538.5 257 539.0 aia ORD 
@. LY: 1.963149 455 562.6 25.8 455 588.4 455 587 — Wie 
N V 1.967393 708 752.0 61.8 708 813.8 708 810 = 3.5 
O VI 1.970911 | 1016969 126 1 017 095 1 017 093 — 3 
F VII 1 973817 1 380 143 232 1 380 375 1 380 421 + 46 
Ne VIII 1 976235 1 798 230 393 1 798 623 = ae 
Na IX 1.978270 2 271 213 626 2 271 839 2 271 832* ra ed 
Mg X 1.980001 | 2799 070 950 2 800 020 2798 834* —186 
Al XI 1.981490 3 381 794 1361 3 383 155 3 383 080 = RD 
Si XII 1.982781 4 019 379 1957 4 021 336 4 021 170* == 160) 
ie) S610 1.983912 4711 780 2689 4714 469 4714 450* == Nts) 


* These terms based on extrapolated values of n* for 3d°D»/,; not fixed by 


Ritz series. 
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Table V 
Calculated and extrapolated values of 2 p?Ps,—3d?Ds), 


Li I sequence. 


v cale. v obs. (Obs.—ceale. ) 

i em + em * em? 
Li J 16 379.14 16 379.12 = 0:02 
Be II 66 121.7 66 117.8 eos 
Beire 147 679.6 147 678.6 —- 10 
CEL; 260 301.5 260 299.6 alo 
N V 403 706.2 403 704.7 i) s S 
ORVE 577 774 577 769 =~ 1) 
avd 782 444 782 487 +43 
Ne VIII = = 
Na IX 1 283 493 1 283 516 + 23 
Mg X 1579 863 1579 904 +4) 
Al XI 1 906 776 1 906 720 — 56 
Si XII 2 264 322 2 264 240 — 82 
P xT 2 652 399 2 652 380 ==) 


where f{(z) is a polynomial. It turns out that such a function fits remarkably 
well for the non-penetrating 2p and 3d orbits in the lithium I sequence. In 
fitting such functions one must remember that the spin orbit interaction term 
must first be subtracted from the observed term value. In the two series 
considered the hydrogenic constants suffice for this purpose. The lower members 
of the sequence can then be utilized to determine the several constants in the 
extrapolation function. Results are shown in Tables III and IV. In all cal- 
culations in this paper the value of R appropriate to the individual atom has 
been used. 

It should be noticed here that only four significant figures are necessary 
in the extrapolation for the 3d terms and five for the 2p terms, due to the 
fact that the extrapolation is run on A. This fact alone makes such methods 
of extrapolation attractive. A further check with experiment is possible if we 
compare the calculated and observed values of the 2 p?Ps, —3d7Ds), transi- 
tion. The results are given in Table V. 

The fact that this method of extrapolation reproduces the minima in the 
curves of n* vs. z and at the same time leads to correct values for the 
2p*Ps,—3d"Ds, transition leads one to believe that some reliance can be 
placed on it. It is possible to fit power series containing terms in z and 1/z 
to these data. Usually one finds, however, that either the differences, term 
obs.—term calc., depart more and more from small numbers or else the slope 
of the curve obtained by plotting 


2p Py, —3 "Ds, — (T2» — Ts) 


against z (where 72, and Tq represent the hydrogenic term values) shows 
constantly increasing deviations. In fact from the nature of such series one 
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might expect such deviations to occur for too wide an extrapolation unless 
the form of the extrapolation formula had an actual theoretical background. 

In fitting the two polynomials given in Tables III and IV to the data, the 
data for the ions through Na IX was considered to be essentially correct and 
tables of third or fourth differences in z?/ A constructed. The higher differ- 
ences are quite sensitive to small experimental error. By trial one can then 
construct a final table with the constant third or fourth difference which gives 
the best fit for the higher members of the sequence. One quickly finds the 
particular ions which seem to fit the most poorly and hence the ones which 
presumably contain the greatest experimental errors. In this connection cer- 
tain of the standard interpolation formulae, such as that due to NEWTON, can 
be most helpful in determining the final numerical constants. 

One can also derive extrapolation functions for penetrating orbits based on 
the older semiquantum mechanical methods: A is, however, much larger than 
in the non-penetrating case. The result is that one is required to retain so 
many significant figures in the calculation that the derivation (which leads to 
quite different extrapolation functions) appears to be of little practical 
interest. 
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